First results for the Coulomb gauge integrals using NDIM 
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The Coulomb gauge has at least two advantadges over other gauge choices in that bound states 
between quarks and studies of confinement are easier to understand in this gauge. However, pertur- 
bative calculations, namely Feynman loop integrations are not well-defined (there are the so-called 
energy integrals) even within the context of dimensional regularization. Leibbrandt and Williams 
proposed a possible cure to such a problem by splitting the space-time dimension into D — u) + p, 
i.e., introducing a specific one parameter p to regulate the energy integrals. The aim of our work 
is to apply negative dimensional integration method (NDIM) to the Coulomb gauge integrals using 
the recipe of split-dimension parameters and present complete results - finite and divergent parts - 
to the one and two-loop level for arbitrary exponents of propagators and dimension. 
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gauges. 
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I. INTRODUCTION 

Perturbative approach in quantum field theory (QFT) was responsible for several breakthrough ideas in Physics 
and Mathematics. One of them is dimensional regularization 0], i.e., analytic continuation of space-time dimension D 
into an extended domain that allows for complex values. Feynman loop integrals gained a solid theoretical foundation 
and renormalization process became simpler than it was (where one had to use cut-offs and so on). Of course, this is 
only a partial picture of it all, say, the covariant side of the coin. 

In algebraic non-covariant gauges J2|, on the other hand, dimensional regularization was able to control divergences, 
e.g., in the light-cone gauge, but the results were not physically acceptable. In other words, double-poles did appear 
in one-loop integral calculations and Wilson loops did not have the correct behavior p|. These problems were first 
overcome with the advent of what is known as the Mandelstam-Leibbrandt (ML) prescription Q . More recently, we 
have shown that in the NDIM approach we do not need to invoke any kind of prescription to perform Feynman loop 
integrals in this gauge ||. 

Among the non-covariant gauges we also have the Coulomb gauge (often referred to as the radiation gauge) where 
confinament Q and bound states |Q are easier to deal with, ghost propagator has no pole and unitarity is mani- 
fest. However, in such a gauge, no further insight has beeen achieved with the standard dimensional regularization 
technique, because it presents a gauge boson propagator of the form 
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which generates loop integrals like, 



d D q 



q 2 {q-p) 2 ' 



(1) 



(2) 



where bold face letters stand for three momentum vectors. 

The integral over the fourth-component (in Euclidean space) or zeroth-component (in Minkowski's), the so-called 
energy integral, is not defined even within the context of dimensional regularization. Doust and Taylor j8| discussed 
Coulomb gauge loop integrals and presented a possible remedy for this problem in terms of a interpolating gauge 
(between Feynman and Coulomb, see also Q). Leibbrandt |10| , again, and Williams, presented another approach 
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for this ill-defined integrals, a procedure called split dimensional regularization. Both parts, namely energy and 3- 
momentum sectors, need to be separately dimensionaly regularized, that is, one parameter only, D, is not sufficient to 
render the integrals well-defined. To overcome this problem, they introduced another regulating parameter, i.e., split 
the dimensinality of space-time into two distinct sectors, namely, D = A — 2e = lu + p and the divergences contained 
in energy integrals are expressed as poles in p besides the usual ones in terms of u>, that is to say, the integration 
measure is written down as dP q = d w q d p q^. 

In a series of papers [[fl| , Lcibbrandt studied Coulomb gauge integrals to one and two-loop level (with Heinrich) 
and presented results for divergent parts of several of them. Our aim, in this work is two-fold: show that NDIM is 
the most versatile technique to carry out loop integrals, whether they come from covariant or non-covariant gauges; 
and to present complete results for the Coulomb gauge integrals to one and two-loop level for arbitrary exponents of 
propagators and dimension. 

The outline for our paper is as follows: in section II we consider scalar and tensorial Coulomb gauge integrals at 
one-loop level, while section III is devoted to two-loop integrals and in section IV we present our conclusions. In the 
appendix we discuss some technical issues. 



II. ONE-LOOP COULOMB GAUGE INTEGRALS 



To show how NDIM can handle Coulomb gauge integrals with ease we consider in this section one-loop integrals. 
Recall that negative dimensional integration is equivalent to positive dimensional integration over Grassmannian 
variables |f 2| — a property demonstrated by Dunne and Halliday — and for this very reason, propagators are raised 
to positive powers (they appear in the numerator of integrands) and usual variables become Grassmannian ones. 
Another important point is that in the NDIM context it is as simple to work with arbitrary exponents of propagators 
as if we choose particular values for them. This is why we consider the general case. It is also worth remembering 
that for some types of diagrams, e.g., box integrals |T|,[l4|], there are divergences that are not related to space-time 
dimension D, but they appear as poles for particular values of exponents, say i and j, of propagators, yielding 
singularities expressed as, for instance T(i — j), see also [ ]15fl . So, within the NDIM approach we can also trace back 
the origin of divergences. 

The first two integrals we choose to work with are scalar ones, 
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92{i,j,k)=l d D q(q 2 y(q + p) 2 i(q 2 ) k , 
where generating functions for these referred Coulomb gauge integrals are, 

Gi= J d D q cxp [-aq 2 - 0(q + p) 2 - 7 q 2 ] , 
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with D = oj + p = 4: — e and d D q — cPq d p q4, in Euclidean space, following the split dimension recipe of Leibbrandt 
et al. 

Completing the square we can easily carry the integration out to get, 
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where Ai = a + (3 + 7. 

Taylor expanding both expressions in (w and (m) we get the NDIM solutions for gi and #2 by solving systems of 
linear algebraic equations. 



The system of linear algebraic equations for the first integral is given by a 5 x 8 matrix, 

X13 + Y u = i 

-X"l23 + Y 25 = j 

X 2 + Y 3 = k 

Y123 = —Xyi — w/2 

Y 45 = -X 3 - p/2 

with five equations and eight "unknowns" , corresponding to the various summation indices coming from the Taylor 
and multinomial expansions. They are solvable only within the lower quadratic 5x5 dimension matrices. There 
are a grand total of 56 possible square matrices of this type (i.e., 5x5) from which 36 yield relevant non vanishing 
and workable solutions while the remaining 20 yield a set of trivial solutions (i.e., the related systems do not have a 
solution). We know from our previous works (see for instance |13(| ) that all those non- trivial solutions will generate 
power series of hypergeometric type, known as hypergeometric functions |lq ]. Moreover, all of them are related by 
analytic continuation, either directly or indirectly. In our present case, namely integral gi, there are triple as well as 
double series, among which we choose to consider only the simplest ones, 

A[AC] r . ,,_ f A,[AC] ^ /p 2 \" 123 {~l) n ?(-i\n 123 ){k + u;/2\n 12 )(D/2 + j\n 23 )(l-i-p/2\n 123 ) 
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where the superscript "[AC]" means analytic continuation (to positive dimensional region) and we define the shorthand 
notation tiab = riA + n B, while (x\y) = (x) y — T(x + y)/T(x) is the Pochhammer symbol and 

/i 4 ' [AC1 = * D/2 {~pl) l {v 2 y +k+D/2 {-]\] + k + io/2)(-k\j + k + D/2){-i + o + w /2|t -a-u/2- D/2 - j - k), 

(10) 

where o = i + j + k + D/2. Observe that the above result is valid for negative j, k. Among the 36 possible series this 
is the only one that has the form S(- • .) a + b + c ; where the • • • stands for the specific kinematical configuration. 
For the case of double series we have, e.g., 

B, [AC]/. ■L.\_ f B,[AC\ v 2 ^ (P 2 Y 12 (-<7\n 12 )(-k\n 2 )(u;/2 + k\ ni )(l - k - a - D/2\n 2 ) 
where 

f B. [AC] = ^/2 (p 2 r( _ 2|ff)( _ j|(j)(cj/2|A;)(A: + a + D/2| _ 2a _ k _ D/2)i (12) 

and the result is valid when i,j are negative. 

This hypergeometric series representation is 4-fold degenerate. Here we mention an important point in the process 
of analytic continuation to positive dimension and negative values of exponents of propagators referred to above. The 
result for the negative dimensional space region for g^ is in fact given by a sum of two terms - the second one being 
also 4-fold degenerate - however, when we perform the analytic continuation this second term vanishes because it 
contains a factor of the form (forgetting about the minus sign) 
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which always vanishes (see also |ll|) since p ^ by definition and where, as usual in NDIM approach we make use of 
the Pochhammer's symbol property (a\b) = ( — l) fc /(l — a\ — b). 

To close this part of the computation, we just mention that of course there are other hypergeometric series which 
represent the same Feynman integral in other kinematical regions, e.g., there is a double series of the form, 
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that is, one of the series (with summation index b) has unit argument and can be recast as a 3 F 2 (...\l). 
These are just the few different manners in which we may write down the result for the integral g\ . 



The second integral is easier than the first, and its result is also degenerate, i.e., there are a total of five 4x4 
systems to be solved of which one has no solution and the remaining four, after properly summed give the same 
result, yielding 
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which after analytically continuing to positive D becomes, 
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Next we consider Feynman integrals in the Coulomb gauge with tensorial structures. Again, as we treated in jlTJ] 
the case of covariant gauge we show here how NDIM can handle these Coulomb gauge tensorial integrals in a similar 
manner. Let, 



and 
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g 4 (t,j, k,m) = J d D q (^(q + p) 2j (q 2 ) fc (2q • p)"\ 
so that, after some algebraic manipulations, we eventually get the result, 
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where the set of parameters {3} = {03, 63, C3; e3, fa} for the hypergeometric function 3-F2({3}) is given in the table. 

We must observe here that for ( |l5| ) the exponent k > always, and in ( |l6| ) the exponent m > always, and these 
must not be analytically continued into the region of negative values, whereas the exponents i,j do follow the usual 
analytic continuation process to get the final result for the integrals. 

From our previous work jlTJ on NDIM approach to tensorial integrals, we know that the best solution for such kind 
of integrals is a truncated hypergeometric function because it contains all the cases of interest in the same formula: 
scalar, vector and arbitrary tensor rank. The hypergeometric function above is clearly truncated for positive integers 
k. This result, among the five possible hypergeometric series representations of such integral, is the only one that is 
a truncated series for even and odd values of propagator exponent k, since it assumes only positive values. 

Finally, the result for the tensorial integral with three propagators, 

# C1 ( J ,^A ; ,m) = 7r B / 2 (-2) m (p 2 r'M|-p/2)(a'+c,/2| J -a0(-jV0H-fc-p/2|-J-^/2) 3 ^ 2 ({4}|l), (18) 



where o J = a + m = i+j + k + m + D/2. Note that the result ([Tq) contains the previous one, (O) in the particular 
case when k = 0; it is valid also when i,j, k are negative and m positive. The five parameters {4} = {04, 64, C4; e^, fi} 
are given in the table, and clearly the hypergeometric function is also truncated for even and odd positive integers m. 
The well-known hypergeometric function 3^2 is defined by the series, 
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so, when we refer to 3-F2 (--• |1) we are meaning the above series. For more details about hypergeometric functions the 
reader is referred to, e.g., |l6[ ]. 



Parameters 


3^ 2 ({3}|1) 


3^2({4}|1) 


a 


-k/2 


-m/2 


b 


1/2 -k/2 


1/2 -m/2 


c 


3 + w/2 


j + lo/2 


e 


1 + i + j + D/2 


l + i+j + k + D/2 


f 


I- i-k- D/2 


1-i-k-m- D/2 



Parameters for hypergeometric functions in equations (O) and (|l§|). 



III. TWO-LOOP COULOMB GAUGE INTEGRALS 

As far as we know, NDIM is the only approach where Feynman integrals in different gauges, covariant and non- 
covariant alike, can be neatly performed, without reference to any special prescription to handle peculiar non-covariant 
singularities in the boson propagator. In the usual covariant gauges several calculations were carried out, e.g., one-loop 
n-point function [18[ , scalar integrals for photon-photon scattering in QED fiM and genuine Jl9| two-loop three-point 
integrals. On the non-covariant side, we have gotten an important original result: Light-cone integrals in the NDIM 
context do not need the famous ML-prescription [Q to circumvent the gauge dependent singularities S, as well as 
avoiding other features which turn the calculation cumbersome - such as using partial fractioning J20| (a mandatory 
feature there) and integration over components. 

Coulomb gauge two-loop integrals can be treated within NDIM methodology as well. Consider, for example, 

I 1 (t,j,k,m) = Jd D qd D r (q 2 Y(r 2 y (p - r - <z) 2fc (r 2 )"\ (19) 

h(i,j,k,m) = J d D qd D r (g 2 r(9-0 2j (r 2 ) fe (p-q) 2m , (20) 
which can be generated by, 

ll= I d D q d D r cxp [-aq 2 - fir 2 - y(p - r - qf - Or 2 ] , (21) 

T 2 = f d D q d D r exp [-aq 2 - [3(q - r) 2 - 7 r 2 - 0(p - q) 2 ] . (22) 

Following the usual steps of NDIM |£1[ , we get for the first integral a 6 x 11 matrix for the system of linear algebraic 
equations to be solved 

-^123 + ^12467 = i 
X\3 + Y 2 368 = j 
-X*123 + ^13578 = k / 9 o\ 

x 2 + r 45 = m ' [Z6) 

X\2 + ^12345 = — W/2 

^3 + Y 678 = -p/2 

which generates 462 (6 x 6) possible hypergeometric series representations for the integral in question. Of these, 216 
have solutions in terms of hypergeometric series whose variable is either z — p 2 /p 2 or z~ x . Among these power series 
the simplest ones are double series which we consider in more detail. Of course, there are also series of the following 
form, 

- z g+» iz -i r +er(...) ~ ^ (^ 1 ) 6+c+e r(-) 

^ a\b\ del TL..Y ^ a\ b\c\e\ T(...Y K ' 

a,b,c,e=0 K ' a,b,c,e=Q K ' 

which can only be convergent if z — z^ 1 = 1. Since this is a particular case, where p 2 = p 2 , we will not study it. 
Let us consider the solution written in terms of double hypergeometric series, 

7^( U)k ,m) = .°( p y' P r £ (^> 12 )(-H^)(- + -/ 2 I^)(;- m ~f'-° /2 M (K) ni \ (25) 

nn=o ni\n 2 \{l - i - k - m - D/2\n2){w/2\n 12 ) \PtJ 

where a" = a' + D/2 = i + j + k + m + D and P A is a product of Pochhammer symbols, 

pjAC] = j-_ .| . + k + D / 2 )(-k\i + k + D/2)(u/2\m)(i + k + D\m){-j\j - o")(o" + m + D/2\j - a"), (26) 

where the exponents of propagators i,j,k must assume negative values. There is another double series, in the other 
kinematical region, where |p 2 /p 2 | < 1, namely, 
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h ^3,k,m)--K {p} P B ^ n 1 \n 2 \(l-a"-u J /2\n 2 -n 1 ){u J /2\n 1 ) {& ) ' 
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where, 

p[AC] = ^^.^^i _ m _ w / 2 )(-fc|j + fc + mH-.D/2)(u;/2|m)(i + ifc + .D| - k -D/2)(p/2\k + w/2). (28) 

The second integral is much simpler than the former, in that the hypergeometric scries representations involved are 
all summable. Summing them is an easy task and the result can be written in terms of gamma functions, 

4 AC] (i,j, k, m) = n D (p 2 ) a " (-i\ - p/2)(a" + w/2| - 2a" - u)/2)(-k\2k + w/2)(-m|2m + w/2)(-j\i + 2j + k + D) (29) 
x(-i-j -k-D/2-p/2\i + p/2)(j + k + D-p/2\ -k-u/2). 

This is a 36-fold degenerate result, i.e., of the overall 56 (5 x 5) systems, 36 of them have non-trivial solutions, 
which after summed (see e.g. pl[ ) and analytically continued give (|2^). It is important to note that this result allows 
negative values for (i,j, k,m) and positive ones for D. For negative m see appendix. 

IV. CONCLUSION 

Splitting the space-time dimension D in the dimensional regularization context into energy-sector and momentum 
sector, each with a specific regularizing parameter, it was possible to Leibbrandt et al to perform perturbative 
calculations in the Coulomb gauge at one and two-loop level. However, the calculations are very involved and they 
were able to present explicit results only for the divergent parts of the integrals. On the other hand, using NDIM we 
showed here that we can calculate complete results for the same integrals, and not only that, they did not have to 
be carried out separately, In our approach we can consider several of them at the same time, because we leave the 
exponents of propagators arbitrary, the integrals being either scalar or tensorial. 
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APPENDIX A: SPECIAL CASES: EXTRACTING POLES 

To make things a little more illuminating we consider in this appendix some technical issues relevant to the results 
for sample special cases. 



1. One-loop. 

Let us consider for instance the particular case where the exponents of propagators in the integral (Hm are i = j 
—l,k = —2, to = 2. The result for this integral is obtained from eq.(|l§|), 
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observing that now a' = D/2 — 2. We proceed as usual in dimensional regularization, taking D = 4 — 2e, u> = 3 — e, 
p = 1 — e and Taylor expanding around e = 0, to get 
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where 7^ is the Euler's constant. 

Now we turn to the theory of hypergeometric functions Kg]. When a numerator parameter is a negative integer 
the series is a truncated one. This is exactly our case, and the hypergeometric function above has only two terms, 
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which substituted into ( |A2j ) yields, 
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We remember that the original integral was contracted with ApiPj, so rewriting 4(p 2 ) e = 4(p 2 ) 2 e 2 
^(PiP l PjP')(p 2 )~ 2 ~ e , we obtain the final result, 
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2. Two-loops. 

Let us consider two particular cases for integral (|2(j): The first one where i = k = —2,j = —l,m= 1 and a second 
one where i = k = — 2, j =m= — 1. 

Observe that in the first case there is one exponent, to, which is positive, so we must not analytic continue it (see 
for instance our previous papers |5|,|l7|). The related Pochhammer symbol (— m|2m + lo/2) was generated by, 
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However it must not be analytically continued since we are interested in the special case where m is positive (to = 1) 
So the result for (p0|), which allows one to take m positive, reads now, 
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Now it is easy to expand (we use the MAPLE V software) around D = 4 — 2e, with a" — D — 4, to obtain the 
poles, double and simple ones, plus hnite part 
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The second special case can be studied using (|29|) directly since in this case all exponents (i = k = —2,j = m = —1) 
are negative, snd (a" = D — 6). Using again MAPLE V software to expand around e = 0, we get a simple pole plus 
finite part, 



4 AC1 (-2, -i, -2, -i) = ^"vr 3 -* (h + l- 2 -¥- 4J f 
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